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Experijnental study resulta for some lattices charac- 
teristics and presented in the paper (The lattices are consis- 
ted of tubular urBnlun( 10-20% enriched) fuel elements 

and light water as a moderator* The experimental study is a 
part of the research program which is underway in the Soriet 
Union on deyelopnent of different type reactors for low power 
nuclear stations /It 2/. 

A special paper section considers syatematization of 
data on critical masses of urani un-water lattices .The seoiempl- 
rical critical equations obtained from such systematisation 
are used to analyse experimental data on the temperature reac- 
tlTity coefficients and excess reactirlty* 

I. CRITICAL EIPERIMEPTS WITH UNDISTURBED COLD LATTICES 
Experiaentsil plant* The majority of the experiments 
the results of which are giren below were performed in a 
special critical facility .The critical facility permits to 
perform exporiaenta in cores with a diameter of 600 mm and a 
height of 800 mm for temperature ranging 13 °C to 230^0. 

The main plant element is a stainless steel high 
pressure vessel with a diameter of 700 mm and a height of 
1800 am designed to operate at pressure of AO atai.At the bot- 
tom of the vessel there is Installed a band electrioeJ. heater 
of 75 kWh oapacity.The upper portion of the Tessel Toluae 
serves as a volume compensator in the experiments with hea- 
ting of the lattices. 
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The experlaentB were made with four types of tubular 
fuel elements, basic parameters of which are presented In 
Table I. 

Critical maaeee* The studied lattices assembled in the 
▼essel of the experimental plant were regular, with no locsd 
perturbations and were of form close to a right circular 
cylinder. The vessel dimensions were such that in the experi* 
ments the lattices were surrounded by an infinite three->sided 
water reflector. The upper reflector was formed by an array 
of the fuel element dry ends .Criticality as rule was achieved 
through vessel filling by water by means of a dosing system. 
Simultaneously with the critical mass measurement the value 
of ^ was measured for each critical height ( jO-the reacti- 
vity, h-the core height corresponding to a given critical 

radius -r ^). The measurements enabled values of the upper 
cr ju 

an«^ bottom water reflector summary savings ( O ) and 
M^/K oo to be determined in the one group diffusion approxi- 
■ation.The experimental data on the critical charge e, dy? /dh 
^ sire given in Fig ■ I and 2 . 

Disadvantage factor for thermal neuti^ns. A disadvanta- 
ge factorC the average thermal neutron flux in a moderator 
to the average flux in a fuel ratio) in the diffusion appro- 
ximation can be written as d=d^ + ^ term due to fuel 

self-shielding. d2 is a term determined by the neutron flux 
non-uniformity in a moderator. 

For fuel elements similar to used d.j«? Ijdg is deteraii- 
ned by a fuel element construction and a lattice pitch, and 
by the thermal neutron diffusion length in water. In complex 
lattices the calculation of dg is not sufficiently reliable, 
ftnH therefore it is required an experimented, verification of 
calculational methods. 

The disadvantage factor was measured in the central 
cell ©f lattices studied from activation of a copper foil 
with a thickness ^ -80 jU- accompanied by the induced 

activity seasursaeat with a system of two end-window / -coun- 
ters .The results of the disadvantage factor measurements at 
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20®C for the fu^l elements of three types are presented in 
Table litas are the calculational data. 

There is a yery small disadvantage factor in the regu- 
lar lattices composed of M-10 and M-20 type fuel elements. 
Howeverflt is assumed the fuel elements of these types to 
use in concrete reactors as fuel assemblies. Therefore water 
gap effect on the disadyantage factor have been studied for 
one of the lattices composed of the type M-10 fuel elements 
!»275)»l?he lattice was assembled with a water gap in 
the centre of it which had a width from 8 to 30 m.Tha disad- 
vantage factor was determined in a plane geometry. An area 
of regular lattice in which the "water gap" effect was still 
observed was taken as the fuel assembly "size". The dependence 
of the disadvantage factor upon the water gap width is given 
in Pig. 3. 

Comparison of experiments with a theory. The critical 
charges and thermal neutron disadvantage factor values obtai- 
ned from experiments were compared with the values calculated. 
The calculations were made on the basis of a spherical harmo- 
nic method, an algorithm of which was described in reference 
/3/*The Immediate calculations were conducted with using of 


proper programs on a digital computer as the 18-group calcula- 
tion in d p, -approximation. Kinetic corrections were obtai- 
ned from a approximation. As it should be expected, the 
experimental and calculated values of are in good agree- 
ment for systems with the ratio ~Ycr j hcV '^0,5. 

Observed discrepancy of the calculations and the experiment e 
at the values ~y^CT / from 0,5 indicates the inade 
quacy of the accepted calculational model for describing of 
the space-energy neutron distribution .Such improvement as an 
account for heterogeneous effects in the resonance area and 
the multigroup thormalization lead only to a slight improve- 

ta.'CJLXtt w.fc vaax? JL 4m.a. v o w xlu ^ ck u Jui/aa 
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tloB With using of appro xlaat ions of oore high order makes 
it permltable a better agreement of the calculational and 
the experimental results In a wider range of the values 

Ycr /^cr 

II.TEMPERATUEE REACTIVITY COEPPICIENTS 
In this paper section results are presented of tempera- 
ture effects study for static conditions (for uniform hea- 
ting of a core and reflectors) tThe measurements were conduc- 
ted in critical assemblies the cores of which were assembled 
from the U-I0>H-20 and K-I7 type fuel elements .The critical 
assemblies composed of the M-IO and M-20 type fuel elements 
were in their characteristics near physically homogeneous 
systems, yet the critical assemblies composed of the K-17 type 
fuel elements were of pronounced heterogeneous structure. 

The h— 10 and U-20 type fuel elements permitted heating to 
100°C only.The K-I7 elements permitted heating to 300®C. 

Disadvantage factor dependence upon temperature. In wa- 
ter-moderated reactors with pronounced heterogeneous core 
structure temperature changes in the disadvantage factor give 
the main contribution to the quantity (C K ^ /(^t and, hence, 
to the teoqperature reactivity ooeff xcient .For several latti- 
ces assembled from the £-17 type fuel elements the disadvan- 
tage factor was studied as a function of the core temperatu- 
re and the ratio J^h/ /^S ( V H^o ^ 

The results of the measurements of the dlaadvantage 
factor temperature dependence for three lattices in the 
range of 20^C -200^0 are presented in Flg.4*The continuous 
lines are the calculatlonal results obtained for the variab- 
les by the method described in section 1 of the paper. 

As it is evident from the curves the calculatlonal data 
are in a qualitative agreement with the experimental results, 
however, a quantitative agreement is inadequate .A dashed line 

Is experimental results processed in a - approximation 
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for a cylindrical block using the experinental teaperature 
dependence of the thernal neutron diffusion length in water 
A/. Such an approxlnation is in better agreement with experi- 
ment and therefore the results were used to calculate dU. /dh 
for the lattices studied. ^ 

Temperat ure effect. In several undisturbed lattices, 
which differed one from another for the most part by the deg- 
ree of their heterogeneity, the effect of temperature on reac- 
of system was studied .The oritioal assemblies were 
heated by an external heat Bourse and thereactivlty change was 
measured as a function of teaperaturs. 

In the physicsQly homogeneous corea(the case of the M-lo 
and M— 20 type fuel elements), the temperature reactivity coeffi- 
cient for which is negative in the whole of temperature range, 
reactivity measurement in heating was performed by a eubcriti- 
cal reactor method. Point of the method is as followsi 

The reactor kinetic equation solution with the external 
neutron source (S) can be expressed in terms of external detec- 
tor counting rate (han^, where N-the external detector counting 
rate, ^ -the counter effectiveness , n-the neutron density 
in a reactor) .The solution for a subcrltical reactor with the 


initial conditions. 


is of the fora 


— 0 cvn/ol/ ■ 

sLie-r 


The quantity L ^ ^ is determined firom the critical 
reactor experiments .Per a critical reactor with the initial 


conditions* 



cj ^ C- 


Mtro 


cLcl 

dv L -S 

Iv-o 
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Heglecting by 1 in compariaon -with 
cing tbe notation 
and (4-) 

fin, = 


6 

I 


/,1 


6ind introdu- 


can be found from (3) 

V - t 

? n . 


2 : 

L 


Thus the equation which ^relates the reactor eubcritica- 
lity to the detector counting rate may be written conclusive- 
ly as ^ 




N (3) 

A neutron source is placed near the critical reactor and a 
curve of counting rate increasing is obtained as a function 
of time. The product « 5^ equals the tangent of the dip 
angle for the straight line corresponding to reactor power 
increasing. 

The present method of measarements requires highly 
accurate organization of experiments. In particular a number 
of detectors uniformly located over critical assembly periphe- 
ry are required, and a bachgroxond must be taken into account 
aepeclally in measurements in case of deep subcriticality. 

In the heterogeneous assemblies which used the K-I7 
type fuel elements at first the reactivity was increasing in 
heating and only after achieving of a certain temperature the 
reactivity was decreasing. Therefore in heating of these 
lattices another method of measurementa was used. In the initi- 
al (taZO^C) the critical assembly was brought exactly to 
criticality. After this the heater was set off and the system 
was heating by ^'^O^Cfthe temperature was stabilizing, and 
the reactor period was measured .Then the system was brought 
to criticality by inserting In the core of absorbers which 
featured little or no effectiveness changes with temperature, 
and the system was heating by ^ts3 +I0^C again. In the experi- 
ments rode placed Inside tubular fuel elements and 


A similar method of the negative reactivity measurement to 
apply for determination of the absorbing rod effectiveness 
was proposed In reference /5/< 

359 
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Btainleas steel tubes enveloping all fuel elements were 
used. This way the measurements were made on the ascending 
branch of the curve p =f(t*^C)* Simultaneously with temperatu- 
re coefficient measurement calibrating of the absorbers was 
going on. The measurements on the descending branch of the 
curve were performed by means of the calibrated absorbers or 
from reactor period measurement as the system was getting 
cold. 

In Pig. 5 the results of the experiments with heating 
of the critical assemblies are presented. 

lo illustrate a leakage contribution to temperature 
effect the curves of heating are shown in Pig. 5 a and 5b for 
the cores assembled from the K-I7 type fuel elements and 
differing one from another with a radius .The dashed lines 
are results of temperature curves calculations .A discussion 
of the calculational and experimental results is given in 
the next section of the paper. 

II I. SYSTEMATIZATION OP EXPER1B4ENTAL WATER 
MODERATED CRITICAL ASSEMBLIES 

Derivation of empirical critical equations. An applica- 
tion of the low group approximations to small water modera- 
ted reactors (k > I»I5j B > 5*10 cm ) meets with a number 

Oo 

of difficulties; 

a) To determine the moderation kernel fourie form It is 
required to account not only for the second moments but mo- 
ments of high order an experimental measurement accuracy of 
which is inadequate. 

b) The extrapolation lengths of small reactors are com- 
parable to core dimensions while the methods of their experi- 
meuta]!. determination have large errors. 

Prom above there was proposed a method by Lanteov M.N., 
one of the u thors of the paper t to obtain semiempirical 

critical equa- 
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tl«ns through systejiati.aUon of experijaental water moderated 
critical asaemblles as follows: 

1) As a characteristic of core multiplicating properti- 
es a calculatlonal value of (k -1) «as used. 

2) As a geometrical factorywhich determines the neutron 
leakage dependence upon core dimensions ,tJtie core surface to 
core volume ratio (s/v) was chosen. 

3) Through systematiaation the dependence ® »f(k -1) 
was found with k^ determined as the secondary \eutron 
number per an absorption .Por all critical assemblies a common 
scheme of calculations was used* 

The systematization with respect to the perameter s/v 
was performed on the base of the present paper data and 
experimental reeulta published before on uranlum~water 
homogeneous and heterogeneous critical assemblies of various 
geometry with /6-9/. Critical equations were obtai- 

ned for various sets of critical assemblies* 

The systematization of heterogeneous critical assemb- 
lies includes only ones for which experimental data could be 
obtained on disadvantage factors and a fraction of U^^^fis- 
sions • 

CyllndricsQ critical assemblies were preliminarily ana- 
lysed ^ith respect to nanrying-out of the condition s/v = 
const as the core dimensions ratio )is changed •Critics! 

assemblies whicbi satisfied the condition within an error of 
no more than - were jncluded in the systematization.lt is 
turned out that the critical assemblies with 0,4 ^ r/h / 1 
satisfied the coiidition«The heterogeneous assembly radius 
was determined from the summary area of core cells. 

In Pig. 6 the results of the systematization of critical 
water moderated assemblies are presented graphically. The ana- 
lytics! resets of the systematization reduce to the follo- 
wing empirical critical equations. 

I) For spberical homogeneous critical assemblies with 
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no reflector 

I = 3/r -0,0M+0,I66( k„-l) 

2)Por apberioal water reflected hoBogeaeoue critical 
assemblies 


I = 0,038 + 0,239 (k ^ - 1 ) ( 5 ) 

»«•> 0,23 4 (k„-l) i 0,92 

3)For cylindrical homogeneoue critical aseemblies with 
no reflector 

I » 2/r + 2/h * 0,053 -f 0,185 -1) (6) 

^)Por cyliodrical homogeneous critical assemblies with 
a full water reflector 

I =.O. 0 « + 0,26 (k„-l) ( 7 ) 

5)Por cylindrical nomogeneoua critical asoembliee with 
a three-sided reflector 


with 


^ =0,04 + 0,246 (k^- 1 ) 


0,15 < (k 


0,93 


with 


6 )Por cylindrical heterogeneous critical assemblies with 

V H a,0 _ 

and a tbrse-sided water reflector 
f =0,026 + 0,243 (k - 1 ) ( 9 ) 


Thus in the considered range of (k^-1) the all criti- 
cal equations is linear and may be written in the following 
general form: 

f-O + COj+Oj) (k^-1) (10) 

where 

C,j- a neutron core leakage characteristic, is a function 
of and core geometry; 

Cg- an independent neutron core leakaee: 

0^-a generalized reflector characteristic independent upon 
both core properties and core geometry. 

For the large values of (k^ -l)the factor s is determined 

for the most part by the coefficients C„ anf 5 *p 
359 2 “4 O3. .3 0^ 4, 
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The Cj universality makes convenient a comparison of 
veiriouB reflectors *In case of such universality for other 


reflector materials 


C:^Ha,0 


S'Hao) 

S' u.n 


(II ) 


where v -the critical size of an unreflected critical 
0686/01)171 -the differences in critical sizes of a 

critical asBembly unreflected and with c given reflector. 

To obtain the calculational dependence s «f( k -1) 
the equation 1-p 

p (12) 

was solved with a given £ and a variable value of Xo 
exrapolation length. To calculate the value of the non-leakage 
probability ( |o ) at slowing-down process the experimental 
slowing-down kernels in water /I0-15/ were used, the Pouri 
form of the slowing -down kernel being determined from 
experimental moments. 

The neutron non-leakage probability calculated from 
the data of reference /lO/ is as 

p = om3+ -fT- atctcuT 8,^se> -^6,5Er +3,56 -ICTlV-e.S- 10 B ^ ^ 

+ 1,9 10 ^ 6 ® 

if the data of references /II-I5/ afe used 

P =-050 -F0 3o 8,18 & - 'IS.S R?' + 3.S6 -10^ dO 

' ^ H9I0= B»] (U) 

In Pig»7 there are given the results of comparison of the 
aistejaatization data with the calculational predictions of 
the neutron non-leakage probability ( p ) through the use of 
the experimental slowing-down kernel Pourie-fora in water. 

The contributions of thermal neutrons and slowing-down 
neutrons to leokage are separated. 

Prom the figure it is seen that the experimental criti- 
cal assemblies systematization results and the calculation 
of p performed u^lng relations (IJ) and (I4) are in good 
qualitative agreement .In the considered range of ( -1) 

the calculation also predicts the linearity of the dependence 
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6 

V *f ^ ^'^oo calculational and expe- 

rimental dependences ^ = | ( k^-1) coincide if the alowing- 
down kernel described in reference /lO/ amd =con8t=3 cm 

are used. 

Application of the empirical critical equations. The 
empirical critical equations obtained from the syatematizati- 
on for heterogeneous critical assemblies were verified by 
comparing with experimental data on reactivity margins for 
critical assemblies of a volume well over critical. 

Such a reactivity margin may be determined experimen- 
tally by graphical integrating of the experimental curve 
dp/dh, = f where h^^^ changes through the critical assemb 

ly radius vafiation. 

Applied to such an experiment the empirical critical 
equations results in 

^ ( 15 ) 


Ak = 


Cj^+ Cj 


Cr 


where ^ - the minimum value of the cylinder critical radi- 
cr 

us with the uranium elements length given; r-a current radius 
value. 

This relation predicts the quantitative results for the 
reactivity margins of the critical assemblies fore-going 
in the paper with an accuracy of the correction to account 
for the delayed neutron efficiency .Here it is important to 
emphasize that the linearity of the dependence between A K 
and ( ), which follows from the systematizati- 

on results, in moat cases fulfils better than in the limits of 
the experimental errors. This moment is illustrated in Pig.8» 
for one of the critioed assemblies studied. 

In interpreting of the reactivity temperature changes 
in critical water moderated assemblies through the empirical 
critical equations it was considered the followijag physical 

picture: 

^S9 
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a) Tile core neutron leakage changee with temperature 
only due to the water density effect on the slowing ~down 
kernel aosents* 


b) The "Infinite" water reflector decreaeee the tenpe— 
rature increase of the core leakage as well as it decrea* 
ses the total leakage .In this case (y » const). 

_ d kpo -y Ca + C3 

pi jr ^3] 


where ^(i) -the water density. 

(T‘^»-":3c0lTC't,)-y(t)] 


Accordingly 

-t 


[Ci]I'(to)+(^J ['“ 1 , ”^^3] 

In the paper the dependence (i'’c.) was found by 

differentiating the following relation for ^ 


(17) 


'00 


'^oo-vL^’e/+ (1-^)0^^] (18) 

where ^ -the number of the secondary neutrons per fis- 
sloni 

^ -the resonance-escape probability in all materials of the 
core 0^ I Qjji^ -fission neutron absorption fractions in 
the thermal group and in the epithermal group relatively* 

In this case the disadvantage factor temperature de- 
pendence chosen was experimental and in the calculation of 
all effects were taken into account contributions 
of which were less than a contribution of the thermal neutrou 
distribution flattening over a cell due to temperature no mo - 
re than by 10 times* 

In Pig*5 the results of the csd.culation of p - | (i'’c) 
using relations (16) and (17) along with the experimental 
dependences p= f which are with no correction for 

the delayed neutron efficiency are presented. Prom the figure 
it is seen that matching of the experimental and calcula- 
tional curves P = | can be achieved by multiplying the 
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esqperifflental data by approxinately constant factor, wJiich 
ia of the same order of magnitude as the correction for the 
delayed neutron efficiency and equals 1,2 7 1,3* 


The critical experimients reported in the paper were 
performed as an eiperimental basis for the reactor design 
of the small transportable nuclear power limit .Obtained 
03 cperlmental information allowed to improye cslculational 
metliods for such reactors and to lay down a program of 
subsequent physical experiments .The performed systematizati- 
on of the critical assemblies reported to date promoted to 
this in large part. 
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PUEL ELEMENTS CHARACTERISTICS 


Table I 




Units 


Type of fuel eleaentn 


M-IO 

M -20 

k-I 7 

k -20 

I. 

Fuel 


UOg 

UOg 

UO 3 

UOg 

2 * 

Enricbaent 

% 

10 

20 

17 

20 

3. 

1 

U^^^amount per I cm 

s 

0.192 

0,264 

3 , 4 a 

3,10 

f-» 

VJI 

of length 






1 

4. 

Number of rxnga 


one 

one 

two 

three 

5. 

Ring size in a can 

mm 

13,4x2,6 

13,4x2,6 

41,8x2,9 

41.8x2,6 






29,2x2,9 

31,6x2,6 







22,2x2,6 

6. 

C€ui material 


Stalnleea 

Stainless 

Stainless 

Stainless 




steel 

steel 

steel 

steel 
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Table II 


CKJ 


CJl 

vO 

DISAOTAMTAGE 

PACTOR POR THERMAL 

NEUTROHS 



'Sffe of fuel 

oleneata 

I 

M-IO 

M>20 


K-I7 


f ft/f, 

^ ezper« 

^ calcul. 

275 

1,05 - 0,02 

260 

1,04 - 0,02 

1,00 

185 

1,87 *0,08 

1,72 

217 

1,98*0,09 

1,87 

252 

1,92 i 0,09 

2,03 
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Fig. 7* The comparison of oalculational and experimental values of 
1,2-.P- calculation using (13) at and 3om, respectively, 
calculation using (14) at7, *2 and 3cm, respectively. 

5 — contribution of thermal neutron leakage. 

6 — results of systematization of bare cylinders. 



Fig. 8. The dependence of reactivity excess on the ratio change of tho sides 
core surface to its voliune. K - 17 fuel element^ f 
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